Bunching and anti-bunching in electronic transport 
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In quantum optics the (^'-function is a standard tool to investigate photon emission statistics. 
We define a ^"'-function for electronic transport and use it to investigate the bunching and anti- 
bunching of electron currents. Importantly, we show that super-Poissonian electron statistics do not 
necessarily imply electron bunching, and that sub-Poissonian statistics do not imply anti-bunching. 
We discuss the information contained in <?' 2 '(t) for several typical examples of transport through 
nano-structures such as few-level quantum dots. 

PACS numbers: 73.63.Kv, 73.50.Td, 73.23.Hk 



Current noise has long-since been established as an im- 
portant tool for studying the physics of transport through 
mesoscopic and nano-scale conductors^— . The charac- 
ter of the noise is typically assessed by considering the 
Fano factor, the ratio of the zero-frequency noise to the 
current^, and comparing with a Poisson process for which 
the Fano factor is equal to one. Systems with F < 1 are 
described as sub-Poissonian (non-interacting systems fall 
in this class^) and systems which have F > 1 are called 
super-Poissonian. A common interpretation of this com- 
parison is that a super-Poissonian Fano factor indicates a 
bunching of the current's constituent electrons, whereas 
sub-Poissonian values indicates anti-bunching (Fig. [T]). 

In this paper we directly investigate bunching and anti- 
bunching in electronic transport as a phenomenon in the 
time domain through the introduction of a second-order 
correlation function g^(r), analogous to that used in 
quantum optics^—. Within a quantum master equation 
(QME) framework in the appropriate limit, the g^- 
function is seen to be proportional to the conditional 
probability that, given an electron is emitted into the 
collector at time t = 0, a further such jump is observed 
a time t later. Following quantum optics, we identify 



g {2 >{0) > 5 (2) (t) bunching 
<?^(0) < g( 2 \r) anti-bunching 



(1) 



since bunching means that particles are more likely to 
be emitted together than apart, and conversely for anti- 
bunching. By relating our g^-function to the correlation 
function between the current at two different times, we 
clarify the relationship between the g( 2 > -function, (anti-) 
bunching and the Fano factor. 

We then investigate bunching and anti-bunching 
in several widely-discussed transport models in the 
Coulomb blockade (CB) regime (see Fig. [2]). This 
analysis shows that the simple picture relating super- 
Poissonian Fano factors to bunching and sub-Poissonian 
ones to anti-bunching is often an oversimplification, and 
can even be outright wrong. In particular we discuss a 
simple quantum-dot (QD) model which has a Fano fac- 
tor less than one, and is thus sub-Poissonian, and yet has 
g^iO) > 3 (2) ( r ) f° r all r > such that, according to 
Eq. (JTJ) , the electron- flow is completely bunched. We 
also give a model for which the converse is true, i.e. we 
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FIG. 1. Sketch of a) anti-bunched, b) Poissonian and c) 
bunched photon (or electron) emission events. 



find a super-Poissonian Fano factor in conjunction with 
electron anti-bunching. These results mirror the work of 
Singh^ and Zou and Mandet^, who have made similar 
points for quantum-optical systems. 

This paper proceeds by first reviewing the situation 
in quantum optics. We then define by analogy our g^- 
function for transport systems and examine its proper- 
ties, including its relationship with current fluctuations 
and the zero-frequency Fano factor. We then consider 
our concrete examples and conclude. 



I. THE SECOND-ORDER CORRELATION 
FUNCTION IN QUANTUM OPTICS 

The second-order degree of coherence of the electric 
field at position r and times t and t + r; r > is defined 



(r, t)E^ (r, t + r)g(+> (r, t + t)£(+> (r, t)) 
(£(") (r, *)£(+) (r, t)) (#(-> (r, t + t)EM ( r , t + r)) 



(2) 



Here, linear polarisation in direction e is assumed such 
that E^ = e ■ . In the stationary limit and sup- 



2 



pressing the position dependence, we have 

(2) <£?(-)£;(-)(r)£;(+)(r)£;(+)) 

with E (± *> = E( ± )(r,t= 0). 

The next step is a re-formulation of the traditional 
Heisenberg operator definition Eq. ([3]) in the (somewhat 
more flexible) Schrodinger (master equation) picture. 
For example, in the theory of resonance fluorescence^, 
one can express the electric field in the far-field limit in 
terms of the atomic (two-level) operators a± (t) via 



E< + >(r,t)=C(r-ro)eo-_(t-<o), 



(4) 



where C(r) is some position-dependent constant, e the 
appropriate polarisation vector and to = r — tq\/c the 
time-of-flight from atom to detector—. With this substi- 
tution, the stationary g^-function of the field can be 
written purely in terms of atomic degrees of freedom as 



9 {2) (r) 



(<t + <t + (t)<t-(t)<t-) 



<7 + <7_ 



(5) 



where all the constants C cancel. The main technical step 
is now a formulation of the dissipative quantum dynamics 
in terms of a quantum master equation with Liouvillian 
W, p = Wp, for the state p of the atom in Born-Markov 
approximation. The quantum regression theorem then 
allows one to writaii 



Tr{ 



<r_e 



Wt 



C-AstatC+W} 



TrjCT-PstatC+r 



(6) 



where p sta t = p{t — > oo) is the stationary state of 
the atom at large times t. One then introduces the 
downwards-jump super-operator J n through its action 
on an arbitrary density matrix p 

J 7 p = -f(T-pcr + , (7) 

where 7 is the spontaneous emission rate, leading to 



9 {2, (r) 



_Tr{j^e WT J yPBta .t} 
Tr{J 7 p sta t} 



(8) 



This equation makes it clear that, at least from a res- 
onance fluorescence perspective, the g( 2 > function mea- 
sures the correlation between two system jumps sepa- 
rated by a time r. 



II. A SECOND-ORDER CORRELATION 
FUNCTION FOR QUANTUM TRANSPORT 

We consider here transport systems in the infinite-bias 
limit such that the time-evolution of the system density 
matrix can be described by a Markovian master equation 
p = Wp*. We assume that we are only interested in the 
current in a single lead, assumed to be a collector, and 



decompose the Liouvillian as W = Wo + J where, sim- 
ilarly to Eq. ([7j), the jump super-operator J describes 
incoherent transitions of the system. Here J describes 
the emission of an electron into the lead in question. The 
remaining part of the Liouvillian, Wo, describes the evo- 
lution of the system without such jumps (jumps to and 
from other leads are included in Wo). When we discuss 
our examples in the next section, we will give the kernels 
in the full-counting-statistics form W(x) = Wo + Je %x 
with x the counting field such that W = W(x = 0) and 
J=-^W( X )|*=o. 

We then define the (/ 2 )-function for transport in anal- 
ogy with Eq. |5|. replacing the optical jump superoper- 
ators with their transport counterparts: 



9 {2 \r) EE 



Wr)J)) 
((J)) 2 



(9) 



where we have explicitly included the master equation 



propagator f2(r) 



and 



denotes the stationary 



expectation value, ((A)) = Tr {Ap sta ,t} with Wp sta t = 0, 
Pstat assumed unique. The interpretation of (r) is 
analogous to the quantum optical case, and in particular 
the designation of bunching and anti-bunching follows 
Eq. Q). 

Whilst we have defined our transport g( 2 )-function here 
by analogy between master-equation formalisms (i.e. in 
terms of jump operators), result Eq. ((HJ) can also be de- 
rived from a microscopic Heisenberg-picture expression 
similar to Eq. @. This derivation is discussed in the 
appendix. 



A. Properties 

We use a Liouville space ket |0}} = p s t a t and bra ((0| 
such that the trace over the stationary state is written 
((A)) = ((0\A\0)). It is then useful to define the projector 
onto the stationary state V = j0>> <<0| via WP = VW = 
0, along with its compliment Q = 1 — V. This allows 
us to decompose the propagator as f2(r) = V + 1Z(t) 
with irreducible part K(t) = Qtt(r)Q. Using ((JVJ)) = 
((0\JVJ\0)) = ((J)) 2 , we find 



g {2 \r) = l + 



((JK(t)J)) 
((JF 



T>0. 



(10) 



In this long-time limit, all non-zero eigenvalues of W have 
negative real-parts and 1Z(t) decays with r — > 00 such 
that £ (2) (t 00) = 1. 

In the limit of strong Coulomb blockade such that at 
most one excess electron occupies the system at any time, 
we have J 2 = 0, since the exit of an electron from the 
system leaves it empty and a second jump can not take 
place immediately. In this case 



ff ( 2 )(r = 0) = 0. 



(11) 



This leads to the conclusion that, according to the defini- 
tion Eq. (JTJ) , a system in the strong CB regime is always 



anti-bunched, since g^'ir = 0) < g^ 2 \r) as g^'ir) is 
always positive. 

The current noise at finite frequency is defined as 
S(ui) = \j dte luJt ({SI{t) 7 SI(0)}). In the QME approach 
with unidirectional transport, this can be written in 
terms of system quantities as 5 - 



S(u) = ((J)) + ((J [ft(-iw) + K{ilo)} J)}, 



(12) 



with Laplace-transformed irreducible propagator IZ(z) = 
[z-QWQP 1 . The mean current reads {I) = ((J)}. 
Transformation back into the time-domain yield a 12 > 13 



S(t) = (I)6(t) + ((JK(\t\)J)); Vr. 
Thus providing we set r > 0, we can equate 
S(t) 



9 {2) (r) 



1; t>0. 



(13) 
(14) 



Or, in the frequency domain, we have 

r+oo 

S(u) = / dre lu,t S{T) 



= (/) + 2{I) 2 / dre tulT 
Jo+ 



(g {2) (r)-l). (15) 



These results show that, in the infinite-bias limit, the 
<7( 2 )-function can be directly related to the current- 
current correlation function either in frequency or time 
domain. 

Defining the Fano factor as the ratio of zero-frequency 
noise to the mean current, we have 

This result for the Fano factor coincides with the opti- 
cal Mandel Q factor Q — F(Q) — 1 in the limit of long 
counting times T — > oor^ 

We also note that the </ 2 )-function defined here bears 
some resemblance to the waiting-time distributio n 15 ^ 6 

iJe w » T J)) 



+ 2(1) ^ dr (V 2 )(t)-i)(.16) 



w(t) 



((J)) 



T > 0, 



(17) 



but there are two differences: (i) the normalisation is 
different, and (ii) the waiting-time propagator is defined 
with Liouvillian Wo which excludes additional jumps 
within the interval r. The g( 2 )-function, in contrast, con- 
tains the full propagator W. 



III. EXAMPLES 

We now discuss the features of the transport g^- 
function for several master-equation models, widely used 
to describe the electronic transport through CB systems 
such as quantum dots (QDs) or molecules. In each case, 
we consider unidirectional transport in a two lead set-up, 
where electrons tunnel into the system from the left lead 
with the rate T^, and out of the system into the right 
with rate F#. 
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(a) Dynamical channel blockade model 



(b) Double quantum dot 




(c) Anderson model for j = 2 



(d) Anderson model for j — 1 



FIG. 2. Sketch of four of the five transport models discussed 
here: (a) The dynamical channel blockade model; (b) the 
double quantum dot; (c) the Anderson model (j = 2); and 
(d) the Anderson model with only one singly-occupied level 
participating in transport (j = 1). In each case, electrons 
tunnel with the rate Yl (Xr) into (out of) the quantum dot 
systems except when the outgoing rate is modified by the 
factor x. 



A. Single resonant level 

Our first example, the single resonant level (SRL), de- 
scribes the transport through a single level (in a QD, for 
example) in the strong CB regime. Written in a basis 
such that the relevant part of the density matrix is col- 
lected into the vector p = {po, Pi}, where and 1 denote 
an empty or occupied level, the Liouvillian reads 



W SRL (x) 



-T L Tr?* 



(18) 



- rRTL and 



The steady state current is given by (/srl) = r , r 
zero frequency Fano factor is^£ 



Fsrl(O) 



(19) 



Since the rates Tl and are (positive) real numbers 
the Fano factor is always sub-Poissonian, F(Q) < 1. The 

(2) 

corresponding ggj{ L (r)-function is given by 



-( r «+ r ^ and g^(0) = 0. (20) 

(2) 

Starting by zero at r = the ggj^ L -function increases 
monotonously to one (not plotted), which indicates 
strictly anti-bunched electrons. Since, here, the Fano- 
factor is sub-Poissonian and the electron flow anti- 
bunched, the SRL model reflects the naive interpreta- 
tion that sub-Poissonian statistics corresponds to anti- 
bunching. Our second example illustrates that this is 
not necessarily the case, however. 



4 




FIG. 3. Main panel: The second-order correlation func- 
tion 5dc B (t) — 1 for the Dynamical Channel Blockade model 
shown for three values of parameter x = 1,1/3,1/5 corre- 
sponding to a sub-Poissonian (dashed), Poissonian (dotted) 
and super-Poissonian (line) Fano factor. Inset: The (shifted) 
Fano factor Fdcb(O) — 1 for the same model as a function of 

(2) 

x. The triangles mark the points for which the <7dcb( t ) ^ s 
shown in the main panel. Other parameters were: Tl — Tr. 



B. Dynamical Channel Blockade 



The dynamical channel blockade (DCB) modet^ is 
an important model of how interaction effects can give 
rise to super-Poissonian statistics. In its simplest ver- 
sion, the model consists of an empty state and a spin- 
up and a spin-down level, see Fig. [Ua). In the basis 
p = {po, P\., p^}, the Liouvillian reads 



-2T L T R e %x xT R e lx 

w uuu ( x ) = | r L -t r o 
r L o -xT R 



(21) 



The dimensionless factor x changes the rate for outgo- 
ing spin-down electrons relative to spin-up ones. The 
steady-state current of the DCB model is (Idcb) = 
and the zero-frequency Fano factor isi£ 



x(r R +r L )+r L ■ 



DCB 



(0) 



r 



r|(3 



2x + 3x 2 ) 



(xT R + xT L +T L )i 



(22) 



For x< 2TL+r t^r^P ' L+L ' R (x<± in the limit T L =T R ), 
the Fano factor of this model is super-Poissonian, see 
Fig- EH (Inset), corresponding to DCB. 

Our (7j 3 Q B (r)-function for this model reads. 



where 7=N /4r2 + (x - l) 2 T 2 R and 7L = T L (x 2 + 1). 

Fig. [3] shows (s , d cb( t ) — 1) plotted for x=l, 1/3, and 
1/5 for which the corresponding Fano factors are 
F DCB (0)=0.55, F dcb (0)=1, and F DCB (0)=1.41, respec- 
tively. 

The sub-Poissonian case (represented by x=l here) 

(2) 

shows a monotonously increasing <? D c B (r)-function with 
interpretation identical to that of the SRL. In the x = 1/5 
case, the Fano factor is super-Poissonian and the func- 

(2) 

tion <?j3c B (t) — 1 starts off negative for r = 0, but rises 
rapidly through zero and remains positive for large r. 

(2) f2l 

For all times, 5dq B (t) > 5dc B (0), and the electron flow 
is anti-bunched. Thus, at least from the perspective of 
definition Eq. ([1]), we find contradiction in associating a 
super-Poissonian Fano factor with anti-bunching, as the 
flow here is anti-bunched independent of the Fano factor. 
Nevertheless, there is a clear difference between super- 
and sub-Poissonian cases, since the -function of the 
former exhibits a clear maximum and a negative slope for 
large times. The area of the curve then integrates 
to a positive number and the Fano factor is greater than 
one. 

This discrepancy can be arises because there are two 
competing processes occuring. At short times the elec- 
trons must naturally be anti-bunched due to the strong 
CB. At longer times, the DCB effect means however, that 
after one electron has tunnelled through the device, there 
is an increased probability that further electrons will tun- 
nel once the system has had a chance to "recharge", and 
the electrons are indeed bunched, but only over time- 
scales larger than is required for the system recharge. 
This analysis mirrors that of Kiefilich et alM- who stud- 
ied the current-current correlation function for this model 
and obtained a similar result to Fig. [3l 

The dotted-curve in Fig. [3] shows the crossover case 
(here at x = 1/3) where the Fano factor is exactly 
-Pdcb(O) = 1. Here the g D( j B (r)-function profile is qual- 
itatively similar to the super-Poissonian case, but the 
"negative-area" exactly cancels the "positive area" result- 
ing in -Fdcb(O) = 1. From the point-of-view of the g( 2 >- 
correlation function nothing particularly special occurs 
at a value x = 1/3, although from the Fano factor, this 
may appear to be the case. Note that, although at this 
point -Fdcb(O) = 1, the complete statistics of the charge 
transfer are clearly not Poissonian since 3 D c B (t) 7^ 1 is 
not equal to one at all times r. This distinction can also 
be shown by calculating the higher-order Fano factors^, 
which for this model deviate from the Poissonian value 
of one. 



C. Double quantum dot 



9 { uUr) =1 



1 

2x7 



(7L+:r7)e^ r(2rL+xIWR+7) 



- {lL ~x~f)e 2 



-ir(+2r Jj +xr H +r H -7) 



. (23) 



Our third example is the double quantum dot (DQD) 
in the strong CB, which allows us to study the influence 
of internal coherence on the (/^-function. A sketch of 
the DQD is shown in Fig.^b) and the system can either 
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FIG. 4. Main panel: The 5dqd (t~ ) function for the double 
quantum dot (DQD) model with zero detuning such that the 
Fano factor is sub-Poissonian, Fdqd(O) « 0.859 (dashed line); 
and detuning e/Tl = 1/2 such that the Fano factor is super- 
Poissonian, F D qd(0) w 1.174 (solid line). Inset: The DQD 
Fano factor as a function of the detuning e.The triangles mark 
the points for which the ffrxjoM - ! ^ s shown in the main 
panel. Other parameters: T C /T L = 1/4, r R /F L = 1/10. 



be empty or have a single electron in either of the left 
\L) or right \R) states. The DQD system Hamiltonian is 
then H DQD = f (\L)(L\ - \R)(R\) + T C (\L)(R\ - \R){L\) 
with detuning e and coupling Tc- The left dot is cou- 
pled incoherently to the source lead and the right dot 
to the collector. With the system density matrix rep- 
resented as p = {poo, Pll, Prr, Plr, Prl}, the system 
Liouvilliar*22r— reads 



w DQD (x) 
' r L 






e lx T R 













iT c 


-iT c 





-Tn 


-iT c 


iT c 


iTc 


-iTc 


-is - T R /2 





—iTc 


iT c 





is - T R /2/ 



(24) 



The stationary current and the Fano factor of this model 

„„»22 



(^dqd) 
-Fdqd 



4T C T L T R 



16 ' 



±T C T R 



T L (8T C 



r|)' 



4T| 



T c 



(l2e 2 



T 2 
1 C 



(4T L e 2 



R ^ 

As 2 



4T 2 (2r z 



r«)) : 



(4T Z 



^r 2 



4T 2 (2Ti 



r^)) : 



.(25) 



Once the right tunnel rate is much smaller than the 
left tunnel rate the Fano factor can become super- 
Poissonian for finite detunings, see Fig.SJInset). For zero 
detuning the Fano factor is sub-Poissonian, but as detun- 
ing increases, the Fano factor becomes super-Poissonian 



around e/Tl ~ 0.15 and reaches a maximum around 
e/Ti w 0.4, before decaying towards unity for large de- 
tunings. 

In Fig. S] we plot the 5dq D (t) function for two sets of 
parameters: one that gives a sub-Poissonian Fano fac- 
tor (dashed line) and one that gives a super-Poissonian 
one (continuous line). In both cases, coherent tunneling 
between the two dots imprints damped oscillations onto 

(2) 

•9dqd( t ) a ^ finite t, similar to as is found in the waiting 
times w(t) for a double quantum dot^, or in resonance 
fluorescence in quantum optics. This is in contrast to the 
other systems without quantum coherences discussed in 
this paper, where g^(r) either monotonously increases 
towards one or has a single peak and then decays to- 
wards one for t — )■ oo. We therefore expect oscillations 
in </ 2 '(t) to be an interesting tool for an experimental 
detection of quantum coherences. 

Note again, that, as for all systems in the ultra-strong 
CB regime, 3dq D (0) = and so the electron flow is anti- 
bunched, according to definition Eq. (JlJ. Once again 
though, this simple assignment misses the full complexity 
of the situation. Both sub- and super-Poissonian traces 
are qualitatively very similar, suggesting once again that 
the value of the Fano factor is itself, not particularly 
diagnostic of the transport mechanisms occurring in this 
model. 



D. Doubly-occupied Quantum dots 

Away from the strong CB limit, the system can be 
occupied by more than just a single electron at a given 
time. This opens the path for the (/^-function at r = 
to be non-zero, and the electron flow bunched. 

As example of this class of model we consider a system 
with one empty state, j £ {1,2} single-occupied states 
(spin orbitals) and one double-occupied states within the 
transport window (and thus contributing to the current). 
In the basis p = {p Q , pi = J2j Pi(j) 1P2}, the kernel reads 



W 2el (x) = 




I* 

r L - t_r xjT R e lx 

T L -xjT R e lx 



(26) 



where x is a factor which modifies the rate of electrons 
tunneling out of the double occupied state. For j = 2 
this model corresponds to the Anderson Model where 
both spin-up- and spin-down-levels are filled (emptied) 
with the same rate Tl (xT r ). In contrast, for j = 1, 
one of the singly-occupied levels lies above the transport 
window. This can be realised, e.g., with a negative charg- 
ing energy^ and a magnetic field. Figures [^c) and ^d) 
shows sketches of these two situations. 

The steady-state current of this model with Fl = T R = 



6 



—3=1, i=l/2 
—3=2, i=l/5 
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-0.10 
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0.0 
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1 2 3 4 5 6 

FIG. 5. Main panel: The (r)-function for the Ander- 
son model with j 6 {1,2} single-electron channels. With a 
single channel (j = 1) and rate parameter x = 1/2 the Fano 
factor is sub-Poissonian (F2 e i(0) = 3/4, see inset), whereas 
the </ 2 '-function clearly indicates bunching (solid line) for all 
times t. With two channels (j = 2) and x = 1/5 the Fano fac- 
tor is super-Poissonian (F2ei(0) = 9/8) but the </ 2 '-function 
indicates complete anti-bunching. This relationship between 
(anti-)bunching and Fano factor is completely opposite to the 
usual interpretation. Inset: Fano factor F(0) — 1 as a func- 
tion of x for j = 1 (line) and j = 2 (dashed). The triangles 

mark the x for which <?dqd( t ) ^ s plotted in the two cases 
where F(Q) 1. Parameters: F L = T R = T = 1. 



(2) 

j = 2 (positive area under <?2ci( r ) — 1" function). This 
is completely opposite to the intuitive understanding of 
the relationship between (anti-)bunching and the Fano 
factor. 



IV. CONCLUSIONS 

The c/( 2 )- function is a well-known tool for investi- 
gating bunching and anti-bunching behaviour of photon 
emission, but can equally well be used to describe the 
statistics of electron emission in transport systems such 
as quantum dots. The unifying frame in both cases is the 
concept of quantum jumps in Markovian quantum mas- 
ter equations that govern the dissipative dynamics. The 
relation Eq. (|16p between the g^- function and the Fano 
factor clarifies that super-Poissonian statistics in elec- 
tron transport do not necessarily correspond to bunching; 
sub-Poissonian statistics do not necessarily imply anti- 
bunching. As our examples show, single electron trans- 
port through nanostructures such as quantum dots offers 
the possibility to test these relations. An interesting fur- 
ther question would be the features of cross correlations 

(2) 

and off-diagonal g kl - functions defined with two different 
jump operators Jk 7^ Ji- 
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r is given by (J 2e i) = Vf +1 . The g\ j (r)-function is 



&r) =1 + ^— f(ai + a 2 )e-M^-7 2 +2)r 



+ (ai - a 2 )e 



-\T(xj+j+ 12 +2)T 



(27) 



with 72 = \/j 2 (x — l) 2 + 4, ai = 272(1 + x — 3jx) and 
a 2 = -2(2(1 + x) +j 2 x{l + x) + j(4x-3 - 3.t 2 )). 



For r = 0, we have 



3X 



1 



Ajx 



(28) 



which is non-zero and becomes g 2c {(fy = 1 f° r x = \ ■ 

In Fig. [5] we plot (5 2 ci( T ) — ■"•) with this value of x for 
j = l (solid line) and j = 2 (dashed line). Although both 
9^ 2 H T )~ 1 curves start at zero for r = 0, the behaviour 
at finite time r is completely different. In the j = 1 

(2) 

model, the <72 C i( T )~function is negative and shows a min- 
imum as a function of time r. The electron flow is there- 

(2) 

fore bunched. In the j = 2 model, however, the 3 2 ei( T )~ 
function is everywhere positive, shows a maximum and 
indicates an anti-bunched electron flow. Looking at the 
the zero- frequency Fano factor (see Fig. [SJ Inset)) we see 
that it is sub-Poissonian for the j = 1 case (negative area 
under (32ci( r ) — -0" function) and super-Poissonian for 
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Appendix A: A microscopic derivation of Eq. (f9|) 

We first specify the Hamiltonian of our transport sys- 
tem as composed of system, lead, and tunnel-coupling 
parts: H = H$ + -f/i ea ds + Hq?. The system Hamiltonian 
we write as Hg = ^2 a E a \a)(a\, where |o) is a many-body 
system state of energy E a . We assume that the leads can 
be described by the non-interacting Hamiltonian 

Pleads = H ^ds = J2Y1 U ak cl k C ak , (AI) 

a a k 

where c a k is the annihilation operator for an electron in 
lead a with quantum numbers k and energy w Q fc. We 
assume that system and the leads are coupled with the 
single-particle tunnel Hamiltonian 



H T = Y v kac\ a d a + H.C., 



(A2) 



fc,Q 



where d a is the annihilation operator of an electron in a 
localised system state, which we assume to be unique for 
each lead, and where Vk a is a tunnel amplitude. 
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To produce a correlation function analogous to Eq. © , 
we first define for lead a the operator 



(A3) 



This definition is analogous to the expansion of electric 
field operator in terms of the normal-mode anni- 

hilation operators^. Here, the choice of coefficients is 
somewhat arbitrary, but 2iV k * a is convenient. We then 
define our (/ 2 )-function for lead-a as 

(2) ,, _ (C { - ] {t)C { - ] {t + T)Ci + \t + T)C { a +) (t)) 



(C { -\t)c£\t)){C { - ] {t + r)C a +) (t + r)> ' 

(A4) 



Heisenberg's equation of motion for the lead annihila- 
tors reads 



Cka{t) = -~iLO ka Ck a (t) - iV ka d a (t). 



(A5) 



Introducing the Laplace transform f(z) = J °° dte zt f(t) 
and solving gives 



Cka(z) = 



1 



{cka(t = 0) - iV ka d a (z)} 



Z + lUJka 

Summing over k and regularizing we obtain 

~ L Z + lUJ ka. ) 

x {c ka (t = 0) - iV ka d a (z)} , (A6) 

where P is the principal part. We choose the initial state 
of lead a such that c ka (t = 0) evaluates to zero in the 
expectation value (see below). Then 



Ci + Hz)=\2nJ2\V ka \ 2 S( 

{ k 



iz - UJ ka ) 
i 1_ 

iz - u) ka 



d a {z) 



= {£ a {iz) + ip a {iz)} d a {z) 



(A7) 



which defines the frequency-dependent rate, T(uj) = 
2nJ2k l^fea 1 2 ~ ^ka), an d principal part p a (ui). If we 
assume that the rate is constant as a function of its argu- 
ment, r a (u;) = T a , and that the principal part vanishes, 
we obtain the result 



C<-+\t)=T a d a (t) 



(A8) 



These assumptions can be justified by assuming 
that lead a starts in the equillibrium state p° q = 
exp {— /3(i/[" ads — p a N a )} in the infinite-bias limit p a — >• 
— oo. This also enforces the requirement c ka (t = 0)p„ q = 
0. 



Placing result Eq. (|A8|) in definition Eq. (|A4[> and tak- 
ing the long-time limit yields 



analogous to Eq. {5J 



,(2)^- (44( T ) d a(T)d Q } st at 



(da rfa)stat 



(A9) 



If we assume that all leads start in this same infinite- 
bias limit, the time-evolution of the system density ma- 
trix can be described by the Markovian master equation^: 



= ~i [H s ,p\ + J2 Ta - \p d a d a - ^d a S a p\ 

r a Idapdlt - Ijpdldg - ^d^dgpl , (A10) 

aeD 



where the first sum is over source leads, the second over 
drains and where r a are rates. Writing this as p = (Wb + 
Sf)p with jump super-operator J p = T a d a pd a , we see 
that definitions Eq. (| A3|) and Eq. (|A4[) indeed reproduce 
Eq. © of the main text. 
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